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STRESS CONCENTRATION AROUND BROKEN FILAMENTS 


IN A FILAMENT-STIFFENED SHEET* 

By W. B. Fichter 
Langley Research Center 

SUMMARY 

Stress distributions around single and double cuts (groups of adjacent broken fila- 
ments) in an idealized composite sheet consisting of parallel tens ion- carrying filaments 
embedded in a shear- carrying matrix are analyzed by an influence-function technique. 
Static-load concentration factors for two equal collinear cuts are obtained and are com- 
pared with some related results for a single cut. Dynamic- load concentration factors, 
corresponding to sudden breaking of filaments, also are obtained for two of the simplest 
double-cut cases. In addition, matrix shear forces are investigated for some single-cut 
and double- cut cases. Finally, loads in broken filaments are calculated for some single- 
cut cases, and their implications in some current failure analyses of composite materials 
are briefly discussed. 

The interaction between two collinear cuts is essentially local and is roughly con- 
fined to separation distances on the order of the cut length. For closely spaced cuts, how- 
ever, the interaction between cuts is seen to be pronounced. For a single cut, an expres- 
sion analogous to the filament load concentration factor is obtained for the maximuni shear 
force in the matrix as a function of the cut length. The results suggest that for some 
combinations of constituent mechanical properties, the composite might be more suscep- 
tible to filament tensile failure at some cut lengths, but more susceptible to matrix shear 
failure at others. 


INTRODUCTION 

Much effort is currently being directed toward the development of lightweight com- 
posite materials typically composed of plastics and metals in which are embedded high- 
strength fibers. An extensive survey of composite materials research is contained in 


*This paper is based on part of a thesis entitled "Stress Concentrations in Filament- 
Stiffened Sheets" submitted in partial fulfillment of the requirements for the degree of 
Doctor of Philosophy in Engineering Mechanics, North Carolina State University, August 
1969. 



reference 1. In studies of the tensile strength of fibrous composites (for example, 
refs. 1 and 2), it has been observed that failure can be initiated by local cuts or breaks 
in the filaments. A theory for predicting composite tensile strength which recognizes the 
statistical nature of filament strength as well as local cut distribution has been developed 
in reference 2. In the statistical approach, however, the stress concentration effects of 
flaw size and flaw distribution are generally ignored or treated rather approximately. 
Therefore, an analysis in which stress concentration effects are taken into account in 
more detail is desirable, and the present paper is directed toward determining some of 
these effects. 

Analytical investigations of the stress fields around some isolated "cutouts" in an 
infinite flat sheet of parallel tension-carrying filaments embedded in a shear-carrying 
matrix are contained in reference 3. Included in these investigations is an analysis which 
yields the filament load concentration factor for a single straight cut across an arbitrary 
number of adjacent filaments (or, equivalently, a group of adjacent broken filaments). 

More recently, this analysis has been extended in reference 4 to problems of load con- 
centration in composites with two-dimensional arrays of parallel filaments. Another 
logical extension of the analysis is to the problem of two cuts in the same filamentary 
sheet, since neighboring flaws in a filamentary composite can be expected to influence 
one another. 

In the present paper an analysis is given of the stresses in a filament sheet weak- 
ened by two collinear cuts of arbitrary length. For the particular case of cuts of equal 
length, filament load concentration factors are computed for various combinations of cut 
length and distance between cuts. In addition, matrix shear forces and filament tensile 
forces are computed for several cases involving single and double cuts. For the case of 
a single cut, an expression analogous to the filament load concentration factor is obtained 
for the maximum matrix shear force as a function of cut length. Finally, the recovery of 
load by broken filaments is investigated for an idealized model and the results, as they 
pertain to statistical strength predictions such as those of reference 2, are briefly 
discussed. 


SYMBOLS 

d filament center-line spacing 

df filament diameter 

EA extensional stiffness of a filament 

E n Weber function 
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G 


matrix shear modulus 


h effective thickness of matrix 

I n modified Bessel function of first kind 

J n Bessel function of first kind 

i,k,m,n,''> . 

> integers 

Q,r,s J 

K filament load concentration factor 

L n modified Struve function 

M mass per unit length associated with a filament 

N n load in nth filament for influence-function solution (L n of ref. 3) 

P n dimensionless load in nth filament, 

p force applied to each filament at infinity 

p„ load in nth filament 

S>max maximum dimensionless matrix shear force 

S n dimensionless matrix shear force per unit length between nth and (n + l)th 

filaments, — J^Qr 
’ P 1/ Gh 

s n matrix shear force per unit length between nth and (n + l)th filaments 

t time 

U n dimensionless displacement of nth filament, 

u n displacement of nth filament 
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displacement of nth filament for influence-function solution 


V n 

x 

r 

T 

e 


coordinate parallel to filaments 
gamma function 

dimensionless coordinate parallel to filaments, 
dimensionless time, t 
variable of integration 




GOVERNING EQUATIONS 


Insofar as is practical, the notation employed in reference 3 is retained in the 
present analysis. The model is shown in figure 1(a), along with the coordinate system 
and some notation. The analytical model is one commonly used in "shear lag" analyses, 
which does not account for local stress variations near the filament-matrix interface. It 
is composed of parallel tension- carrying members (filaments) embedded in a matrix which 
carries only shear. The filaments are separated by a constant distance d and are num- 
bered from -oo to oo. The coordinate along the filaments is x and the displacement of 
the nth filament at location x is u n (x). The force in the nth filament, taken positive in 
tension, is denoted by p n (x) and is given in terms of u n by 

*>„« - EA It <» 


where EA is the extensional stiffness of the filament. The shear force per unit length 
between the nth and (n + l)th filaments is defined here by s n = ~^u n+ j - u n j. Static 
equilibrium of an element of the nth filament then requires 


dx 


p n ( x ) + s n (x) - s n _ 1 (x) = 0 


or, in terms of displacements, 


EA 


d ^ u n Gh 


dx 


# + T( U n + l- 2u n^n-l) = » 


( 2 ) 


In figure 1(a), filaments -2, -1, 2, and 3 are shown cut at x = 0 for illustrative 
purposes. In general, for two collinear cuts through q and s filaments, let 
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-(k + q) £ n S -(k + 1) and r+l§nSr + s, respectively, denote the broken filaments, 
the two cuts being separated by r + k + 1 intact filaments. (See fig. 1(b).) Then the 
boundary conditions at x = 0 are 


u n (0) = 0 (-k ^ n ^ r;-(k + q + 1) ^ n;n ^ r + s + l) 

p n (0) = 0 (-(k + q) S n = -(k + l);r + lSn^r + s) 


(3) 


For x large, the force in each filament approaches the uniform applied force, denoted 
by p; that is, 

P n (±°°) = p (4) 


For convenience, let 


Pn = P P n 


Ujl P ^EAGh U n 


/G h“ 0 
PyEAd Sn 




(5) 


Then the equilibrium equations and boundary conditions become 

d 2 U n 

— # +U n + l- 2U n + U n .i = 0 

dr 

U„(0) = 0 (ni -(k + q + l);-k = n S r;n = r + s + l) ^ 

P n (0) = 0 (-(k + q) ^ n = -(k + l);r + 1 ? n i r + s) > 

P n (±°°) = 1 


The dimensionless forces and displacements are related by 

dU, 


Pn “ d| 


n 


S n = U n+1 - U n 


( 6 ) 


(7) 


(8) 
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SOLUTIONS 


The mixed boundary value problem defined by equations (6), (7), and (8) can be solved 

by use of the influence-function technique employed in reference 3 to solve the single- cut 

dV n 

problem. The influence functions V n (4) and N n (4) = —^(4) are, respectively, the non- 

dimensional displacement and force in the nth filament when the filament sheet is com- 
pletely cut along 4 = 0 and the zeroth filament is displaced a unit amount at 4 = 0 while 
all other filaments are held fixed at 4 = 0. In terms of V n (4) and N n (4) (L n (4) in 
ref. 3), the dimensionless force and displacement in the nth filament are given by 


OO 

P n (£) = l + ^ N n-i^ U i(°> 


l=-00 


OO 

v 


u„(?) = 4+ ^ v l . i ( i )» 1 (0) 


l=-0O 


(9) 




In reference 3, V n (4) has been found to be 

r7r -24 sin I 

V n (4) = ± \ cos ne e z de 

n v J 0 

Application of boundary conditions (7) to equations (9) yields 


-(k+1) 


r+s 


P n (4) = l+ ^ N n _ i (4) Uj(0) + ^ N n _i(4)Ui(0) 

i=- (k+q) 


-(k+1) 


i=r+l 


r+s 


U 


n (4) = 4 + £ V n -i(4) Ui(0) + £ V n-i(4) U i(0) 

i=- (k+q) 


i=r+l 


j 


r+s 


since Uj(0) = 0 for other values of i, and 

-(k+1) 

0 = 1 + \ N n-i(°) Ui(0) + ^ Nn.iWU^O) 

i=-(k+q) i=r+l 


(10) 


(ID 


(-(k + q) ^ n ^ -(k + l);r + l£nSr + sj (12) 
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which expresses the condition of zero load on the ends of the broken filaments. Equa- 
tions (12) constitute a set of q + s linear algebraic equations in the q + s unknowns 
U n (0). Their solution set can be substituted into equations (11) to obtain expressions for 
load and displacement in any filament. 

However, before the loads and displacements can be calculated, the integral 
representation of the influence functions V n (£) and N n (£) must be evaluated. Only 

N n (0) = — — (0) has been evaluated in reference 3 because the computation of load con- 

d£ 

centration factors, which was the main purpose of that investigation, does not require the 
evaluation of the influence functions for nonzero values of ij. Equation (10) can also be 
written as 

V n ($) = - P cos 2n0 e -2 ^ Sin 9 d 9 (13) 

TI J Q 

Integration of this expression gives 


V n (£) = J 2n (2i|) - iE 2n (2i£) 


where i denotes y/^T, J 2 n is the Bessel function of the first kind of order 2n, and 
E 2 n is the Weber function of order 2n. (See ref. 5.) In terms of functions with real 
arguments, 


v n ($) = (-D r 


I 2n (2?) - L 2n (2£) - I £ (-1) 


n-1 r 
k 


( k - D £ 


1\ >.2n-2k-l 


k=0 


t” + r 


(n § 0) (14) 


in which l 2 n is the modified Bessel function of the first kind, L 2n is the modified 
Struve function (see ref. 5), and the finite sum must be taken to be zero for n = 0. A 
similar expression for V n (£) can be derived for n < 0; however, it is obvious from 
equation (13) that V_ n (ij) = V n (^). Differentiation of equation (14) yields 


N n («) = (-D^n+l* 2 ^ + hn-l^ “ 


t2n 

L2 ” +l(25) + L 2n-‘ <2{) %l /2 r( 2 „ + |) 


n-l 


- i I (- 1 ) 1 


k=0 


(2n - 2k - l)r(k + ^ 

r (2„ + 1 - k) 


2n-2k-2 


(n ^ 0) (15) 


in which the finite sum again is to be omitted for n = 0. From equation (15) it is found 
that 
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(16) 


N n (0) = 

7r(4n 2 - l) 

which is in agreement with reference 3. 

With substitution of the appropriate values from equation (16), equations (12), which 
are merely linear algebraic equations, can be solved for the displacements of the ends of 
the broken filaments U n (0). Then the solutions to equations (12), along with the influence 
functions given by equations (14) and (15), can be substituted into equations (11) to obtain 
the load and displacement at any point in the filamentary sheet. However, because each 
cut may traverse any number of filaments and the two cuts may be separated by any 
number of intact filaments, the general case of two collinear cuts is without symmetry. 
Hence, calculations covering reasonably wide variations of the pertinent parameters 
would entail considerable computational effort. If the problem is simplified, however, by 
requiring the cuts to be of equal length, the most essential features of the two-cut prqblem 
are retained and the computational effort is greatly reduced through consideration of the 
resulting symmetries. 

In what follows, then, it is assumed that the two collinear cuts are of equal length. 

In the analysis of this reduced problem, one of two cases arises, depending upon whether 
the number of filaments between the cuts (henceforth called "interior filaments") is even 
or odd. The analysis of these two cases is presented in the appendix, pertinent results 
being drawn from there as appropriate. 

Filament Loads 

In the case of a single cut across n consecutive filaments (beginning with the 
zeroth filament), the load and displacement in the rth filament were found in reference 3 
to be 

n-1 

P r (£) = 1 + £ N r _i(£) U 4 (0) (17) 

i=0 


Ur(S) = « + £ V r _ t (€) Uj(0) (18) 

i=0 

where the U^(0) terms are the solutions to the n simultaneous equations 

n-1 

0=1 + ^ N r _ t (0) 1^(0) (Olrln-l) (19) 

i=0 
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For the case of two equal-length collinear cuts separated by any number of intact fila- 
ments, expressions for the load and displacement in any filament are derived in the 
appendix. 

In reference 3 the load concentration factor for a single cut across n filaments 
was found to be 

n-1 

K n = P n (0) = 1 + Y N n _i(0) Ui(0) (20) 

i=0 

which is the maximum load in the filament just ahead of the cut. For double cuts, the most 
highly stressed filaments are the outermost of the filaments between the cuts and, again, 
their maximum load occurs at | = 0. In the appendix, expressions for the load concen- 
tration factors have been found for cases involving odd and even numbers of interior fila- 
ments. For an odd number 2r + 1 of intact filaments between two cuts, each across n 
filaments, the load concentration factor is 

r+n 

K n,2r+l = 1 + £ j^r+i(0) + N r _ i (0)Ju i (0) (21) 

i=r+l 

where the Uj(0) terms are the solutions to equations (A4) for a specific value of 
for an even number 2r of interior filaments, 

r+n-1 

K n ,2r = l+ £ [N r+i (0) + N r _ i _ 1 (0)Ju i (0) 
i=r 

where the Uj(0) terms are calculated from equations (A10). 

Shear Loads 

The dimensionless shear force per unit length between the nth and (n + l)** 1 filaments 
is given by 

S n (0 = u wl («) - U n (g) (23) 

In either a single-cut or double-cut configuration, the greatest matrix shear force is 
between a broken filament and an unbroken one. In the case of a double cut, the unbroken 
filament is an interior one. 

To find the matrix shear forces in a single-cut problem, the solutions to the boundary 
load conditions (19) are substituted into equation (18) to obtain displacements, and those 
results are used in equation (23). In a double-cut problem, either equations (A4) or (A10) 


n; and 
( 22 ) 
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are solved, depending on the number of interior filaments, and the solutions are substi- 
tuted into the second of either equation (A3) or (A9), respectively, and then equation (23) 
is used. 

Calculations of matrix shear forces have been made for some single-cut and double- 
cut cases. The calculations are restricted in each case to the matrix area between two 
filaments which experiences the greatest shear stress. The solutions are: 

For single cut across one filament: 

S 0 (S) = f ^ - 2[l 0 (2£) - L 0 (2£)J - |[L_ 1 (2^) - I 1 (2^ (24) 

where 

s 0 (0) = - l 

For single cut across two filaments : 


S l© - - L 0<2«] + |(2 + |i)[Vl<2?> - Il<2{)] - 

where 

Si(0) = - 

For single cut across three filaments: 


s 2 U) = - 







20 + 


237 





where 


(25) 


(26) 


S 2 <°) - - ft 

For two cuts, each across one filament, with one interior filament: 

Sod) = §^2 + ij[lo(2{) - L 0 (2«)] ♦ |(l + -4-^1 (27) 

where 



10 



For two cuts, each across two filaments, with one interior filament: 


^ + ?)[v 2 «> - l o< 2 «] * f f - ^)M« - 

-!r-^-i|( 1 + il + ^)&> (25) - L o <2 «3 

+ f( 3 + fl + ff)[ L -i <2 «>- 'i«!l - k- 3$) 

where 


(28) 


c m\ — 55 65 tt 
S ° ( } " 2(6512) 

and where the bars are used to identify double-cut shear forces. A comparison of equa- 
tions (27) and (28) will indicate how rapidly the complexity of double- cut shear force cal- 
culations grows with cut length. 

RESULTS AND DISCUSSION 

The solutions presented in the preceding section have been employed to investigate 
various phenomena associated with the breaking of filaments, such as static and dynamic 
load concentration in filaments, matrix shear force growth with cut length, as well as 
decay of matrix shear force and recovery of load by broken filaments as functions of axial 
distance from the cut. 


Load Concentration Factors 

Tensile loads .- In reference 3 the static load concentration factor for a single cut 
across n consecutive filaments was found to be 


_ 4 • 6 • 8 • • • (2n + 2) 
n " 3 ■ 5 • 7 • • • (2n + 1) 


(n = 1, 2, 3, . . .) (29) 


In the present paper, filament load concentration factors for two collinear cuts K n m 
each across n filaments, separated by m intact filaments, have been calculated for 
n from 1 to 8 and m from 1 to 16. The results are presented in table I and are plotted 
in figure 2, along with the single-cut factor K n for purposes of comparison. It might 
be noted here that K n is the asymptotic value of K njin for large m, that is, for large 
distances between the cuts. It can be seen in figure 2 that the interaction between collin- 
ear cuts is essentially local, being confined to separation distances which are on the order 
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of the cut length. For cuts which are farther apart, the results do not differ conspicuously 
from their asymptotic, single- cut values. Thus, when two cuts are separated by distances 
which are greater than their lengths, each cut can be considered independently with little 
loss of accuracy. 

The values of K n>m for closely spaced cuts, however, increase rapidly when m 
becomes small. Hence, K n)in for m small is of particular interest, since it is asso- 
ciated with states of severe load concentration. In order to examine the relative severity 
of double cuts and comparable single cuts, double- cut factors for closely spaced cuts are 
plotted in figure 3 along with factors for single cuts of comparable total length. In terms 
of the notation employed here, K n j and K n> 2 are compared with single-cut factors 
denoted by K 2 n , K 2 n+ j, and K 2 n+ 2 - These particular values of the single-cut factor 
are employed because they represent either the same total number 2n of broken fila- 
ments or the number 2n + 1 or 2n + 2 which would be broken if the double cuts sepa- 
rated by one or two intact filaments, characterized by K n> i or K n ^, were to coalesce. 

It can be seen in figure 3 that the load concentration factors for the double cut with 
one interior filament K n? i is greater than those for a single cut across an equal number 
of filaments when more than two filaments are broken. Also, the factors for the double 
cut with two interior filaments K n 2 are greater than those for a single cut across an 
equal number of filaments when 15 or more filaments are broken. Thus, it appears that 
the longest cut does not necessarily induce the most severe stress state. 

This result suggests that a composite structural design criterion based on residual 
strength of the composite in the presence of a single cut of prescribed length could be 
unconservative, even though only cuts of lesser length were present. 

In reference 3 dynamic load concentration factors were obtained for cases in which 
a single cut is suddenly introduced in a stretched filament sheet and, in a separate analy- 
sis, an apparent upper limit of 1.27 was found for the dynamic overshoot (ratio of maxi- 
mum dynamic and static load concentration factors). To investigate the possibility of a 
departure from the trend for a single cut, dynamic load concentration factors have been 
computed for two of the simplest double cut cases. The dynamic analysis was carried 
out by inserting the appropriate inertia term in equation (2) and utilizing the Laplace 
transform theoi’y in the manner of reference 3. The results are presented in figure 4, 
where some single cut results are reproduced for purposes of comparison. For the two 
double-cut cases, in which totals of two and four filaments were suddenly broken, one 
remaining intact between them, the dynamic overshoots were found to be 1.22 and 1.23, 
respectively, and are in the range of previously obtained single-cut results. The fact 
that additional double-cut calculations have not been made because of their rapidly 
increasing complexity precludes the drawing of sweeping conclusions concerning the 
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overall behavior of the double-cut dynamic overshoot; however, the present limited cal- 
culations suggest that departure from the range of single- cut values is unlikely. 

Shear loads .- Although the present composite model accounts for matrix shear 
stress in only an approximate manner, it is believed that matrix shear-force calculations 
are of at least qualitative interest. A result which was not obtained in reference 3, but 
which can be extracted from the analysis presented there, is the magnitude of the most 
severe shear force at a single cut as a function of the number of cut filaments, which 
might be thought of as a shear-force concentration factor S max . For a cut which starts 
at the zeroth filament and severs n filaments, the most severe shear force is given by 

^max = |s n (0)| =u 0 (0) (30) 

where S max is defined as the absolute value of the peak shear force, and where Uq(0) 
is the nondimensional displacement of the end of the zeroth filament, which varies with 
the number of cut filaments. Solution of equations (19), (18), and (23) for the first six 
values of n yields the following values of 


n 

1 

2 

3 

4 

5 

6 


Smax ( see also e( l s - (24), (25), and (26)): 


g 

max 

0.785 

1.178 

1.473 

1.718 

1.933 

2.126 


Inspection of these values in fractional form indicated that they conformed to the 
expression 


4 C _ ( 2n - !)• 

ir D max - 9 9 To 

2 2n “ 2 [(n - I)’] 


(31) 


By using equations (27) and (28), a few results were also obtained for the shear load 
concentration factor S max for the two simplest double- cut configurations in which a sin- 
gle intact filament separates the cuts. Solutions for 2 and 4 broken filaments yield 


n 

2 

4 


®max 

0.841 

1.342 
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Unlike tensile load concentration factors for closely spaced cuts, . these limited results 
suggest that the shear load concentration factors for double cuts are less than those for 
single cuts with the same number of cut filaments. 

A further consideration in comparing shear and tensile load concentration factors 
is their relative rates of increase with increasing number of cut filaments. In figure 5, 
the shear concentration factor S max for single cuts is compared with the tensile load 
concentration factor K n for values of n up to 13. The shear force can be seen to 
increase more rapidly than the tensile force, particularly for small values of n. Both 
K n and S max are unbounded as n approaches infinity. However, their relative mag- 
nitudes for large n can be determined by writing equation (29) in the form 

Tjr _ 2 2n n T .(n + I)'. 
n (2n + 1): 

and forming the ratio 

\ s max n(2n+ l)[(2n)i) 2 
Kn = 2 4n - 1 (n + l)(n*.)4 


By the use of asymptotic formulas for the factorial function, it is found that 


q 

°max j 

K n 


(n - oo) 


(33) 


Equations (32) and (33) indicate that the maximum dimensionless matrix shear force and 
filament tensile force are of comparable magnitude for all values of n. They are approx- 
imately equal for large values of n, although, as noted, the maximum shear force initially 
increases more rapidly than the tensile force. This difference in slopes for small n 
suggests that some composite materials could be more susceptible to tensile failure of 
the filaments when only small flaws are present, but might be more susceptible to matrix 
shear failure when weakened by larger flaws. 

The longitudinal variation of the shear loads was also investigated. In figure 6, the 
most severe matrix shear forces for single cuts across one, two, and three filaments (see 
eqs. (24), (25), and (26)) are plotted as functions of the axial distance from the cut. Also 
shown are results for the two simplest double-cut configurations (see eqs. (27) and (28)) 
in which a single intact filament separates the cuts. The curves are normalized with 
respect to their peak values, which can be obtained from the table of amplitudes in fig- 
ure 6 along with the pertinent equation from equations (5). It can be seen that the peak 
values increase with cut length. The normalized curves differ from one another in their 
rates of decay with axial distance from the cut, the decay generally being slower for 
longer cuts. It might also be noted that the double -cut shear loads decay somewhat 
more rapidly than those for the single cuts. 
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Loads in Broken Filaments 


For a single cut across n consecutive filaments (starting with the zeroth fila- 
ment), the load in the ith filament is given by 

n-1 

Pi(£) = 1 + £ Ni_ k (4) U k (0) 
k=0 

where the n broken filaments are identified by i = 0, 1,. . . n - 1, and each U k (0) is 
one of a unique set associated with a specific value of n. The U k (0) values are obtained 
by solution of equations (19). 

Loads in the broken filaments have been calculated for a single cut across one, two, 
and three filaments. The results are plotted in figure 7, where it can be seen that the 
length over which a broken filament carries less than a particular fraction of its load can 
vary considerably with the number of broken filaments. The 90-percent load-recovery 
line has been added for later reference. With slightly more labor, similar results can 
be obtained for various double-cut configurations; however, results for double cuts are 
not expected to exhibit significantly different trends. For this reason, and because illus- 
tration of the effect of cut length on filament load recovery is the main purpose of the cal- 
culations, loads in broken filaments have not been calculated for double cuts. 

Ineffective Length Calculation 

In reference 2, the ineffective length of a filament is defined as that portion of a 
broken filament which supports less than 90 percent of the load applied at infinity. Studies 
cited in references 1 and 2 suggest that the ineffective length, based on a load-recovery 
fraction of 90 percent, of a single broken filament embedded in a matrix can range from 
less than one up to several hundred filament diameters, depending on the geometry and the 
mechanical properties of the constituents. Calculations have been made in reference 2, 
for example, of ineffective length as a function of filament volume fraction Vf, and E/G, 
the ratio of filament Young's modulus to matrix shear modulus. Results are presented in 
the form of a family of curves, each member of which corresponds to a particular filament 
volume fraction. 

In the case of a composite containing a cut, an additional parameter which would be 
expected to influence the ineffective length is the number of broken filaments. With the 
aid of an idealized model, results from the present study can be used to obtain an indica- 
tion of the influence of this parameter. A typical cross section of this model is shown in 
the sketch. For the present ineffective length calculations, it is assumed that the thick- 
ness of the sheet is equal to the filament diameter, and that in the matrix between any two 
adjacent filaments the shear stress is constant at a given axial location. Then for this 
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model the filament volume fraction is given by v f = where df is the diameter of a 

1 4d 

filament and d is the filament spacing. 

In terms of the pertinent mechanical and geometrical parameters, the axial distance 
along a filament is given by 


x = 





which for the idealized model becomes 


X_ _ 7T /X* 

df 4tfGv f 5 


(34) 


In order to obtain estimates of ineffective length changes due to the breaking of additional 
filaments, it is necessary to determine from figure 7 the values of denoted by at 
which the broken filaments of interest have recovered a stipulated fraction of their far- 
field load. These values of £, along with appropriate values of E/G and Vf , are then 
substituted into equation (34) to obtain values of the nondimensional ineffective length, 
denoted by x'^/df. 

This procedure has been carried out for two values of vj and a range of values of 
E/G, the filament being considered ineffective wherever it supports less than 90 percent 
of its far-field load. The cases of cuts across one, two, and three adjacent filaments have 
been considered. The results are presented in figure 8. In the case of the cut across 
three filaments, the 90-percent recovery figure pertains to either of the two outer 
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filaments in the cut rather than to the middle filament which would have yielded even 
greater values of ineffective length. 

The mathematical model employed in the analysis of reference 2 consists of a sin- 
gle filament encased in a thin layer of shear-carrying material (binder) which, in turn, 
is embedded in an infinite body to which are assigned the average stiffness properties of 
the composite material. Hence, interaction between neighboring filaments is ignored. In 
addition, shear stresses in the average material are assumed to decay in a negligible dis- 
tance from its interface with the binder material. Because of these basic differences 
with the present model, the results of reference 2 are not plotted in figure 8. It should 
be noted, however, that the results of reference 2 yield a single curve for each filament 
volume fraction, whereas the present results yield a family of curves, each member of 
which corresponds to a particular number of broken filaments. In addition, the analysis 
of reference 2 predicts a more rapid recovery of load by broken filaments than does the 
present analysis. As can be seen in figure 8, large changes in filament ineffective length 
can result from varying the number of broken filaments. The results suggest that this 
effect should perhaps be accounted for in statistical studies of the strength of filamentary 
composites. 


CONCLUDING REMARKS 

An analysis has been made of the loads and deformations in a filament- stiffened 
sheet which has been weakened by two collinear cuts. Also, some results for the case of 
a single cut, in addition to those previously reported, have been obtained. It has been 
found that significant interaction between collinear cuts is largely restricted to cases in 
which the distance between cuts is no greater than the cut length. It is seen that two 
closely spaced cuts can cause greater stress concentration than a single cut across a 
comparable total number of filaments. Limited calculations of dynamic stress concen- 
tration factors for suddenly introduced collinear cuts support the earlier conclusion that 
dynamic effects are not of great importance in filamentary composites of the type 
investigated. 

Calculations of loads in broken filaments show significant changes in ineffective fila- 
ment length with cut length; thus, statistical strength analyses of composite materials 
should consider the incorporation of a flaw-size parameter. 

Calculations for single cuts of various lengths show that maximum shear forces grow 
more rapidly with cut length than maximum tensile forces. Such calculations might be 
useful in determining whether a composite material is more susceptible to shear failure 
or tensile failure. 
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The present analysis is based on linear small-deflection theory of a single sheet, 
whereas filamentary composites usually are many filaments thick and are subject to var- 
ious nonlinear effects such as plastic deformations, large deflections, and straightening of 
the filaments. Although these effects would be expected to exert some perturbing influ- 
ence, the present results would be expected to remain qualitatively descriptive of load 
concentration effects in filamentary composites of the type investigated. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Langley Station, Hampton, Va., June 30, 1969. 



APPENDIX 


STRESS CONCENTRATION FACTORS FOR TWO 
EQUAL-LENGTH COLLINEAR CUTS 


In this appendix attention is restricted to problems in which the two collinear cuts 
are of equal length. This restriction is effected by setting q equal to s in equa- 
tions (11) and (12); this procedure yields 


-(k+1) 


r+s 


Pn(£) = l+ £ N n _ i (S)U 1 (0)+ £ N n . i (€)U i (0) 

i=(k+s) 


-(k+1) 


i=r+l 


r+s 


(Al) 


u n («) = | + I, V n _ i (?)U i (0)+ £ V n . i («)U i (0) 

i=-(k+s) 


and 


■(k+1) 


i=r+l 


r+s 


0 = 1+ ^ N n _i(0) Ui(0) + Y, N n-i(°) u i<°) 

i=- (k+s) i=r+l 


(-(k + s) § n g - (k + l);r + 1 Sn^r + s) (A2) 

One of two cases arises, depending on whether the number of intact filaments between the 
cuts (interior filaments) is odd or even. The case of an odd number of interior filaments 
is treated first. 


Odd Number of Interior Filaments 

When the two cuts are of equal length, the displacements are symmetric with respect 
to the line which equally divides the group of interior filaments. In this case it is con- 
venient to take the zeroth filament as the line of symmetry so that k = r in equations (Al) 
and (A2) and the total number of interior filaments is 2r + 1. With Uj(0) = U_j(0) 
because of symmetry, equations (Al) and (A2) become after some manipulation 
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and 


r+s "> 

Pn(«) = l+ J [ N n + i(^) + N n _i(^Ui(0) 
i=r+l 

r+s 

U n (l) = €+ J jv n+i (£) + V n _i(^Ui(0) 

i=r+l 


(A3) 


r+s 

0 = 1+ ^ [N n+i (0) + N n _ i (0)Ju i (0) (r+l?nsr + s) (A4) 

i=r+l 

The maximum load will occur adjacent to the cuts (£ = 0) in the two outermost interior 
filaments (n = ±r). Hence, from the first of equations (A3), the stress concentration fac- 
tor for two cuts across s filaments separated by 2r + 1 filaments is 

r+s 

K s,2r+1 = p r(°) = 1 + £ j^r+i(0) + N r _ i (0)]u i (0) (A5) 

i=r+l 

where the Ui(0) terms are calculated from equation (A4) for specified values of r 
and s. 


Even Number of Interior Filaments 


In this case, the interior filaments are assumed to be identified by -kJn^r, 
which with r = k - 1 becomes -klnSk - 1. Equations (Al) and (A2) take the form 


P„(5> 


u n (S) 


-(k+1) 



k+s-1 



1+ I 

N n- 

■ id) Ui(0) 

♦ I 

N n -i(l) 

Ui(0) 

i=-(k+s) 



i=k 



-(k+1) 



k+s-1 




v n - 

id) Ui(o) 

♦ l 

V n -iU) 

Ui(0) 

i=-(k+s) 



i=k 




(A6) 


and 
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-(k+1) k+s-1 

0 = 1+ Y, N n-i(0) Ui(0) + Y N n _i(0)Ui(0) 
i=- (k+s) i=k 

(-(k + s) § n S -(k + l);k ^nSk + s- 1)) (A7) 

For this configuration, the axial line of symmetry is located between the filaments num- 
bered 0 and -1, which means that 

u _i(0) = U^O) (A8) 

After substitution from equation (A8) and some manipulation, equations (A6) and (A7) 
become 


and 


r+s-1 

P n^ = 1+ I §Wl(£) + N n -i(Sj u i(°) 

i=r 

r+s-1 

U n (4) = ^+ Y |7n+i+l^) + v n-i^)]Ui(0) 
i=r 






(A9) 


r+s-1 

0 = 1 + Y [ N n+i+l’(°^ + Nn-i^ 0 )]^ 0 ) (r s n S r + s - 1) (A10) 

i=r 


In this case the maximum load occurs for n = -r, r - 1. Then with n = r - 1 in the 
first of equations (A9), the stress concentration factor for two cuts across s filaments, 
separated by 2r filaments, is 

r+s-1 

K s,2r = p r-l(°) = 1 + Y £W°) + N r _i-i( o g u i(°) (All) 

i=r 

where the Uj(0) terms are calculated from equations (A10) for specified values of r 
and s. 
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TABLE I.- LOAD CONCENTRATION FACTORS FOR TWO 
COLLINEAR CUTS OF EQUAL LENGTH 


Number of 
filaments 
between 
cuts 

1 

2 

3 

4 

5 

6 

7 

8 
9 

10 

11 

12 

13 

14 

15 

16 


Load concentration factors when the number of 
filaments in each cut is - 


1 

2 

3 

4 

5 

6 

7 

8 

1.714 

2.359 

2.964 

3.543 

4.102 

4.646 

5.177 

5.698 

1.412 

1.788 

2.142 

2.481 

2.808 

3.126 

3.436 

3.740 

1.368 

1.690 

1.989 

2.272 

2.543 

2.805 

3.061 

3.310 

1.353 

1.654 

1.928 

2.185 

2.430 

2.666 

2.895 

3.118 

1.346 

1.636 

1.897 

2.140 

2.370 

2.590 

2.803 

3.009 

1.342 

1.626 

1.879 

2.113 

2.333 

2.542 

2.745 

2.940 

1.340 

1.620 

1.867 

2.095 

2.308 

2.510 

2.705 

2.893 

1.338 

1.615 

1.859 

2.082 

2.290 

2.488 

2.677 

2.859 


1.612 

1.854 

2.073 

2.278 

2.471 

2.655 

2.833 



1.849 

2.066 

2.268 

2.458 

2.639 

2.813 


1 

1.846 

2.061 

2.261 

2.448 

2.626 

2.798 




2.057 

2.255 

2.440 

2.616 

2.785 




2.054 

2.250 

2.434 

2.608 

2.775 





2.246 

2.428 

2.601 

2.766 





2.243 

2.424 

2.595 

2.759 






2.420 

2.590 

2.753 
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(a) Double cuts across filaments -2, -1, 2, and 3. 
Figure L- Coordinate and notation systems. 
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(b) Collinear cuts across q and s filaments, separated by r+k+1 intact filaments. 
Figure 1.- Concluded. 




Figure 3.- Comparison of filament load concentration factors for single and double cuts. 
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Figure 1 .- Loads in broken filaments near a single cut. 







